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Abstract 

During 2010 the LHC delivered pp collisions with a centre-of-mass energy of 7 TeV. 
In this paper, the results of comprehensive studies of missing transverse energy as 
measured by the CMS detector are presented. The results cover the measurements of 
the scale and resolution for missing transverse energy and the effects of multiple pp 
interactions within the same bunch crossings on the scale and resolution. Anoma- 
lous measurements of missing transverse energy are studied, and algorithms for 
their identification are described. The performances of several reconstruction algo- 
rithms for calculating missing transverse energy are compared. An algorithm, called 
missing-transverse-energy significance, which estimates the compatibility of the re- 
constructed missing transverse energy with zero, is described, and its performance is 
demonstrated. 
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1 Introduction 

Neutral weakly interacting particles, such as neutrinos, escape from typical collider detectors 
without producing any direct response in the detector elements. The presence of such particles 
must be inferred from the imbalance of total momentum. The vector momentum imbalance in 
the plane perpendicular to the beam direction is particularly useful in pp and pp colliders, and 
is known as missing transverse momentum, here denoted Bt- Its magnitude is called missing 
transverse energy, and is denoted B t- 

Missing transverse energy is one of the most important observables for discriminating leptonic 
decays of W bosons and top quarks from background events which do not contain neutrinos, 
such as multijet and Drell-Yan events. It is also an important variable in searches for new 
weakly interacting, long-lived particles. Many beyond-the-standard-model scenarios, includ- 
ing supersymmetry predict events with large Bj- The reconstruction of Bj is very sensitive to 
particle momentum mismeasurements, particle misidentification, detector malfunctions, parti- 
cles impinging on poorly instrumented regions of the detector, cosmic-ray particles, and beam- 
halo particles, which may result in artificial Bj- 

In this paper, we present studies of £ t as measured using the Compact Muon Solenoid (CMS) 
detector at the Large Hadron Collider (LHC), based on a data sample corresponding to an in- 
tegrated luminosity of 36 pb _1 . In Sectionji] the CMS detector is briefly described. In Section|3j 
particle reconstruction algorithms and identification requirements, together with the basic sam- 
ple selection requirements, are given. In Section ij the different algorithms for evaluating Bj 
are presented. In Section [5J methods for identifying anomalous Bj mismeasurements from 
known detector artifacts are described. In Section [6j the scale and resolution are measured 
using events containing photons or Z bosons. The degradation of the resolution due to the 
presence of additional soft pp collisions in the same crossing as the hard scatter ("pile-up") is 
presented. In Section|7|we present distributions from physics processes containing genuine Bj- 
In Section |HJ an algorithm, called "Bt significance", which is the likelihood that the observed 
B t is due to resolution effects, is described, its performance in jet events is demonstrated, and 
its efficacy for separating events containing a W boson decaying to either an electron and a 
neutrino or a muon and a neutrino from multijet backgrounds is shown. Conclusions are given 
in Section|9j Finally, in the Appendix, the optimization of the parameters used in the correction 
for the detector response is described. 



2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diam- 
eter, providing a field of 3.8 T. Within the field volume are the silicon pixel and strip tracker, 
the crystal electromagnetic calorimeter (ECAL), and the brass/ scintillator hadron calorimeter 
(HCAL). Muons are measured in gas-ionization detectors embedded in the steel return yoke. 
In addition to the barrel and endcap detectors, CMS has extensive forward calorimetry. 

CMS uses a right-handed coordinate system, with the origin at the nominal interaction point, 
the x-axis pointing to the centre of the LHC, the y-axis pointing up (perpendicular to the LHC 
plane), and the z-axis along the anticlockwise-beam direction. The polar angle, 9, is measured 
from the positive z-axis and the azimuthal angle, (p, is measured in the x-y plane relative to the 
x-axis. Transverse quantities, such as "transverse momentum" (pj), refer to the components in 
the x — y plane. The magnitude of pj is pj. Transverse energy, Ej, is defined as E sin 6. 

The electromagnetic calorimeter (ECAL) consists of nearly 76 000 lead tungstate crystals, which 
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3 Data sample selection and particle reconstruction 



provide coverage in pseudorapidity \rj\ < 1.479 in a barrel region (EB) and 1.479 < \rj\ < 3.0 
in two endcap regions (EE). A preshower detector consisting of two planes of silicon sensors 
interleaved with a total of 3Xo of lead is located in front of the EE. The ECAL has an energy 
resolution of better than 0.5 % for unconverted photons with Ej > 100 GeV. 

The HCAL is comprised of four subdetectors, a barrel detector (HB) covering \rj\ < 1.3, two 
endcap detectors (HE) covering 1.3 < \t]\ < 3.0, two forward detectors (HF) covering 2.8 < 
< 5.0, and a detector outside of the solenoid (HO) covering \tj\ < 1.3. The HCAL, when 
combined with the ECAL, measures hadrons with a resolution AE/E 100% \JE [GeV] © 
5%. In the region \r\\ < 1.74, the HCAL cells have widths of 0.087 in pseudorapidity and 
0.087 rad in azimuth. In the (f],(p) plane, and for \rj\ < 1.48, the HCAL cells map onto 5x5 
ECAL crystal arrays to form calorimeter towers projecting radially outwards from close to the 
nominal interaction point. At larger values of \t]\, the size of the towers increases and the 
matching ECAL arrays contain fewer crystals. 

The muons are measured in the pseudorapidity window \rj\ < 2.4, with detection planes made 
of three technologies: Drift Tubes, Cathode Strip Chambers, and Resistive Plate Chambers. A 
global fit of the measurements from the muon system and the central tracker results in a pj 
resolution between 1 and 5%, for pj values up to 1 TeV. 

The inner tracker measures charged particles within the \fj\ < 2.5 pseudorapidity range. It 
consists of 1440 silicon pixel and 15 148 silicon strip detector modules and is located in the 3.8 T 
field of the superconducting solenoid. It provides an impact parameter resolution of ~ 15 f/m 
and a pj resolution of about 1.5% for 100 GeV particles. 

The first level (LI) of the CMS trigger system, composed of custom hardware processors, uses 
information from the calorimeters and muon detectors to select, in less than 1 fis, the most 
interesting events. The High Level Trigger (HLT) processor farm further decreases the event 
rate from around 100 kHz to ~ 300 Hz, before data storage. 

The calibrations used in this analysis were those available at the beginning of the fall of 2010, 
which are not the final calibrations for this data sample. Improvements and updates to the 
calibrations can lead to small improvements in resolution. 

A much more detailed description of CMS can be found elsewhere (T|. 

3 Data sample selection and particle reconstruction 

The data sets used for the studies presented in this paper were collected from March through 
November, 2010, and consist of pp collisions at a centre-of-mass energy \/s = 7 TeV. An in- 
tegrated luminosity with all subdetectors certified as fully functional of 36 pb _1 was avail- 
able. The detailed selection criteria for the individual data samples used for each study are 
given throughout the text. However, all require at least one well-identified primary vertex 
(PV) whose z position is less than 24 cm away from the nominal centre of the detector and 
whose transverse distance from the z-axis is less than 2 cm, ensuring that particles coming 
from collisions are well contained in the CMS detector. 

The samples used for the studies in this paper are defined through selection requirements on 
reconstructed jets, electrons, photons, muons, and b jets. We describe the basic identification 
requirements used for these particles here. 

Jet reconstruction and its performance in CMS are discussed in detail elsewhere [2j. For the 
analyses described in this paper, jets are reconstructed using an anti-kj algorithm [3j with a jet 
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radius parameter R of 0.5. The energy of a jet is corrected, on average, to that which would 
have been obtained if all particles inside the jet cone at the vertex were measured perfectly 
(particle-level). CMS uses three different types of jets. Calorimeter jets (Calo Jets) are clus- 
ters of calorimeter tower energies. Jet-plus-track jets (JPT Jets) achieve improved response by 
supplementing the calorimeter information with tracking information. Tracks are associated 
with Calo Jets if they are within the jet cone at the PV. The measured momentum of these 
tracks is added to the jet. To avoid double-counting energies, the expected response in the 
calorimeter is subtracted from the Calo Jet if the particle is still within the cone when it im- 
pacts the calorimeter. The response and resolution can be further improved using a global 
particle-flow reconstruction. Details on the CMS particle-flow algorithm and performance can 
be found in Q. The particle-flow technique reconstructs a complete, unique list of particles (PF 
particles) in each event using an optimized combination of information from all CMS subdetec- 
tor systems. Reconstructed and identified particles include muons, electrons (with associated 
bremsstrahlung photons), photons (including conversions in the tracker volume), and charged 
and neutral hadrons. Particle-flow jets (PF Jets) are constructed from PF particles. 

Photon candidates are selected from clusters of energy in the ECAL. They are required to be iso- 
lated. The ECAL energy in an annular region in the fj — (p plane with inner radius 0.06 and outer 
radius 0.4, excluding a three-crystal-wide strip along cp from the candidate, is required to be less 
than 4.2 + 0.006 pj GeV, where pj is the transverse momentum in GeV of the photon candidate. 
The sum of the pjs of the tracks in the same region must be less than 2.2 + 0.0025 pj GeV. The 
ratio of the HCAL energy in an annular region with inner radius 0.15 and outer radius 0.4 to 
the ECAL cluster energy is required to be less than 0.05. The shape of the cluster of calorimeter 
energies must be consistent with that of an electromagnetic shower. In addition, to provide 
strong rejection against misidentification of electrons as photons, the cluster must not match 
any track reconstructed in the pixel detector that is consistent with coming from the primary 
vertex. Photon reconstruction and identification are described in detail in 0. 

Electrons are identified using similar criteria. In addition to similar shower shape and isolation 
requirements, the candidate must match well in both (p and tj to a charged track, but be isolated 
from additional tracks. Electron candidates are also required to be in the fiducial portion of the 
calorimeter (\rj\ < 1.4442 or 1.5660 < \t]\ < 2.5). More details are given in [6]. In addition, 
photon-conversion rejection is used in some of the analyses presented in this paper. 

Muon candidates consist of a track in the tracker which can be linked to one reconstructed in 
the muon system. The candidate must be isolated from deposits of energy in the ECAL and 
HCAL that are not consistent with having been deposited by the muon. The sum of the pjs of 
other tracks within an isolation cone centered on the candidate must also be small. 

Several algorithms for the identification of b jets have been developed flZHEl. Two of them are 
used in an analysis described in this paper. The SimpleSecondaryVertex (SSV) tagging 
algorithm exploits the significance of the three-dimensional flight distance between the PV and 
a reconstructed secondary vertex. The Sof tMuonByPt (SMbyPt) tagger uses the transverse 
momentum of the muon with respect to the jet axis to construct a discriminant. 

The collision data are compared to samples of simulated events that were generated either 
using PYTHIA 6 |9|, with a parameter setting referred to as tune Z2, or with MadGraph inter- 
faced with PYTHIA IfTOl . The generated events are passed through the CMS detector simulation 
based on GEANT4 IfTTl . The detector geometry description includes realistic subsystem condi- 



tions such as simulation of nonfunctioning channels. The samples used in Section 6.5 include 



pile-up. The offline event selection for collision data is also applied to simulated events. 
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4 Reconstruction of E T 



4 Reconstruction of 2£ T 

In general, £ j is the negative of the vector sum of the transverse momenta of all final-state par- 
ticles reconstructed in the detector. CMS has developed three distinct algorithms to reconstruct 
£x : ( a ) PF£t, which is calculated using a complete particle-flow technique |4J; (b) Calo£x, 
which is based on calorimeter energies and the calorimeter tower geometry [Tl2| ; and (c) TC £ j, 
which corrects Calo £ j by including tracks reconstructed in the inner tracker after correcting 
for the tracks' expected energy depositions in the calorimeter |Tl3l . 

PF £ x is calculated from the reconstructed PF particles. PF J2 Ej is the associated scalar sum of 
the transverse energies of the PF particles. 

Calo £ j is calculated using the energies contained in calorimeter towers and their direction, 
relative to the centre of the detector, to define pseudo-particles. The sum excludes energy 
deposits below noise thresholds. Since a muon deposits only a few GeV on average in the 
calorimeter, independent of its momentum, the muon pj is included in the Calo £ j calculation 
while the small calorimetric energy deposit associated to the muon track is excluded. Calo J2 Ej 
is the associated scalar sum of the transverse energies of the calorimeter towers and muons. 

TC £ x is based on Calo £ j, but also includes the pjs of tracks that have been reconstructed in 
the inner tracker, while removing the expected calorimetric energy deposit of each track. The 
predicted energy deposition for charged pions is used for all tracks not identified as electrons 
or muons. The calorimetric energy deposit is estimated from simulations of single pions, in 
intervals of pj and rj, and an extrapolation of the track in the CMS magnetic field is used to de- 
termine its expected position. No correction is applied for very high pj tracks (p T > 100 GeV), 
whose energy is already well measured by the calorimeters. For low-px tracks (pj < 2 GeV), 
the measured momentum is taken into account assuming no response from the calorimeter. 

The magnitude of the £ j can be underestimated for a variety of reasons, including the nonlin- 
earity of the response of the calorimeter for neutral and charged hadrons due to its noncompen- 
sating nature, neutrinos from semileptonic decays of particles, minimum energy thresholds in 
the calorimeters, pj thresholds and inefficiencies in the tracker, and, for Calo £ j, charged par- 
ticles that are bent by the strong magnetic field of the CMS solenoid and whose calorimetric 
energies are therefore in a calorimeter cell whose associated angle is very different from the 
angle of the track at the vertex. The displacement of charged particles with small pj due to the 
magnetic field and the calorimeter nonlinearity are the largest of these biases, and thus Calo £ j 
is affected most. A two-step correction has been devised in order to remove the bias in the 
£t scale. The correction procedure relies on the fact that £j can be factorized into contribu- 
tions from jets, isolated high pj photons, isolated high pj electrons, muons, and unclustered 
energies. The contribution due to unclustered energies is the difference between the £x and 
the negative of the vector sum of the pjs of the other objects. Isolated photons, electrons, and 
muons are assumed to require no scale corrections. 

Jets can be corrected to the particle level using the jet energy correction [21 . The "type-I cor- 
rections" for £x use these jet energy scale corrections for all jets that have less than 0.9 of their 
energy in the ECAL and corrected pj > 20 GeV for Calo £ j, and for a user-defined selection of 
jets with pj > 10 GeV for PF £j. These corrections can be up to a factor of two for Calo £x but 
are less than 1.4 for PF £ j [(14) . In order to correct the remaining soft jets below this threshold, 
and energy deposits not clustered in any jet, a second correction can be applied to the unclus- 
tered energy, which is referred to as the "type-II correction". This correction is obtained from 
Z — > ee events, as discussed in the Appendix. 
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In this paper, distributions involving Calo £ j include both type-I and type-II corrections, those 
involving PF £ j include type-I corrections, and those involving TC £ j are uncorrected, as these 
were the corrections that were available at the time the analyses presented in this paper were 
performed and are the versions used most typically in 2010 physics analyses. As discussed 
in the Appendix, type-II corrections have been developed for PF £ x and can be used in future 
analyses. The optimization of both corrections is also discussed in the Appendix. 

5 Large E T due to misreconstruction 

This section describes various instrumental causes of anomalous £t measurements, and the 
methods used to identify, and sometimes to correct, £ j for these effects. We also examine the 
contributions to the tails of the £ j distribution from non-functioning channels, uninstrumented 
regions of the detector, and particles from sources other than pp interactions. 

5.1 Contributions to E T from anomalous signals in the calorimeters 

The CMS ECAL and HCAL occasionally record anomalous signals that correspond to particles 
hitting the transducers. Anomalous signals in HCAL can also be produced by rare random 
discharges of the readout detectors. Some of these effects had already been observed during 
past test beam and cosmic data taking lfTSl . Detailed studies of these effects have continued 
with the 7 TeV data taking, and are documented in detail in fl6| for the HCAL, and in IfTT) for 
the ECAL. For some types of anomalous energies, the number of affected channels is small and 
the event can still be used in physics analysis after the removal of the anomaly. We refer to the 
removal process as "cleaning" the event. If a large number of channels are effected, "filters" 
instead tag the event as not suitable for use in physics analysis. 

Anomalous energy deposits in EB are associated with particles striking the sensors and very oc- 
casionally interacting to produce secondaries that cause large anomalous signals through direct 
ionization of the silicon. Three main types of noise have been identified in HF: scintillation light 
produced in the light guides that carry the light from the quartz fibres to the photomultipliers, 
Cherenkov light in the photomultiplier tube (PMT) windows, and punch-through particles hit- 
ting the PMTs. While the EB, HF scintillation and HF Cherenkov sources typically affect only 
a single channel, signals generated in the HF by particles that exit the back of the calorimeter 
can affect clusters of channels per event. In the HB and HE, electronics noise from the Hybrid 
Photo Diode (HPD) and Readout BoX (RBX) occurs, and can affect from one up to all 72 chan- 
nels in an RBX. This noise is not related to interactions with particles from pp interactions but 
instead occurs at a low rate and at random times, so the overlap with pp interactions is very 
low at the bunch spacings of the 2010 run. 

The basic strategy for the identification and removal of anomalous signals (cleaning) is based 
on information such as unphysical charge sharing between neighbouring channels in rj-cp and/ or 
depth, and timing and pulse shape information. Each of the calorimeters in CMS measures and 
samples signals every twenty-five ns and several samples are saved with the event record. The 
shapes of the pulses for signals that develop from energy deposits in the calorimeters are dif- 
ferent than those from anomalous noise signals. 

Once a "hit" in an HCAL tower or ECAL crystal is determined to be unphysical, we exclude 
it from the reconstruction of higher-level objects such as jets or £j- We thus arrive at a recon- 
struction of jets and £x that is consistently "cleaned" of anomalous detector effects. Studies 
using simulations of a variety of different physics processes indicate that the amount of energy 
due to particles produced in a pp scattering that is removed is negligible. 
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Some features of anomalous signals can be used to identify events contaminated by them most 
effectively after higher level objects such as jets have been reconstructed. Usually we reject 
events containing these types of anomalies using filters instead of trying to clean them, al- 
though some cleaning is available for PFBj- For example, we usually exclude events with 
HPD or RBX noise affecting many channels from our data samples using the requirements de- 
scribed in [15]. We find that this requirement excludes 0.003% of an otherwise good inclusive 
sample of pp interactions (minimum-bias events). 

The Calo £ j distribution from a data sample that was collected on a trigger that requires a coin- 
cidence in the beam pick-up monitors and scintillators in front of the HF calorimeter (minimum 
bias data) are shown before and after removal of the anomalous signals in Fig. [lj demonstrat- 
ing the effect of the cleaning and filters. A comparison with simulation, which does not include 
anomalous energies, shows good agreement. The effect of the cleaning on the other types of 
£ x is similar. 

The minimum bias triggers used to collect the data shown in Fig. [I] were prescaled for most of 
the data-taking period. Triggers that require large amounts of energy in the detector, such as 
£t triggers and single jet triggers, are enriched in events with anomalous energies. Filters for 
cosmic rays, other non-collision-related sources of high £ j, and other types of anomalies have 
been developed in the context of specific searches for new particle production [181. An example 
of a filter for beam-halo muons, which can produce high energy bremsstrahlung photons in the 
detector, is given in the next section. 

Hardware modifications to mitigate one of the largest sources of anomalous energies during 
the 2010 run, scintillation light produced in part of the light guide reflective sleeves in the HF, 
were implemented during the winter 2010 shutdown period. During this period, the material 
that was producing the scintillation light was replaced with Tyvek. These modifications reduce 
the rate of noise events in HF by an order of magnitude. The HF PMTs will be replaced with 
multi-anode PMTs with flat, thinner front glass during 2013/2014 winter shutdown, reducing 
the noise from Cherenkov light and punch-through particles. To reduce the noise observed in 
HB and HE, HPDs will be replaced by Silicon Photo Multipliers (SiPM) that do not produce 
this type of noise. 

5.2 Removal of beam-induced contributions to E T 

Machine-induced backgrounds, especially the production of muons when beam protons suf- 
fer collisions upstream of the detector ("beam halo"), can cause anomalous, large £j. The 
CMS beam-halo event filter uses trigger and reconstruction-level information obtained from 
the Cathode Strip Chambers (CSCs), a subdetector with good reconstruction performance for 
both collision and non-collision muons [19| and can be used to tag events for removal. The 
geometry of the CSCs make it difficult for beam-halo particles, with mostly parallel-to-beam 
trajectories, to traverse the barrel calorimetry without traversing one or both CSC endcaps. 

The filter can operate in either a "loose" or "tight" mode. The former is designed for high tag- 
ging efficiency at the cost of a modest misidentification probability, while the latter tags only 
well-identified halo candidates and has a smaller misidentification probability. The tagging 
efficiencies and misidentification probabilities have been assessed using simulation. For sim- 
ulated beam-halo particles which impact the calorimeters and produce Calo Bt > 15 GeV, the 
loose (tight) filter is roughly 92% (65%) efficient. The per-event mistag probability determined 
from a simulation of inclusive pp interactions (minimum-bias events) for the loose (tight) filter 
is ~10~ 5 (10~ 7 ). The tagging inefficiency is due in part to halo muons which do not traverse 
enough active layers of the CSCs for a well-measured track to be reconstructed and in part to 
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Figure 1: Calo Bj distributions in a minimum bias data sample without (black dots) and with 
(open circles) cleaning and filters, compared to simulation. Overflows are included in the high- 
est bin. 

muons that do not meet the coincidence requirements of the LI beam-halo trigger. Many of the 
mistagged events are from extremely soft and forward muons (i.e. pj < 2 GeV and \tj\ > 1.7), 
from pion decay or from hadron punch through. 

The CSC-based beam-halo filter was applied to events passing muon triggers which had pj 
thresholds of 9, 11, or 15 GeV, depending on the running period. Beam-halo muons, because 
their tracks do not point towards the nominal interaction point in the centre of the detector, 
in general do not fire the triggers for muons from pp interactions. Thus the beam halo muons 
in this sample are overlaid on events triggered otherwise. This sample therefore provides an 
unbiased comparison of ?t in events with and without a beam-halo muon in coincidence. 
Minimum bias events could have been used as well, but, because the minimum bias trigger 
was prescaled, the number of available events was small. The fraction of halo-tagged events 
for each running period is shown versus the average beam intensity, with an uncertainty of ap- 
proximately 10%, in Fig.[2peft). The fraction of tagged events increases with the beam intensity, 
as might be expected. 

Figure fright) shows the PF B j distribution for two trigger streams. The distribution from 
events recorded by collision muon triggers is shown by the dashed curve while that of the 
subset of these events which met the requirements of the tight halo filter is shown by the red 
inverted triangles. As can be seen, the halo muons that overlapped with these events did not 
disproportionately produce events with large PF 1 j, which indicates that the probability that 
a halo muon produces large Bj in events taken from triggers that are uncorrelated with Bj 
is small. However, events from a trigger on Calo B t with a minimum trigger threshold of 
100 GeV (solid curve), show a substantial fraction identified as halo (blue triangles), since the 
trigger preferentially selects events in which the beam-halo muon has deposited large amounts 
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of energy in the calorimeter. A beam-halo filter is therefore necessary for analyses that make 
use of samples based on this trigger. 

5.3 Contributions of non-instrumented or non-functioning detector regions 

Particles traversing poorly instrumented regions of the detector can be a cause of apparent Bj. 
While generally hermetic, the CMS calorimeter does have uninstrumented areas (cracks) at the 
boundary between the barrel and endcap sections, and between the endcap and the forward 
calorimeters. The gap between the barrel and endcap sections is about 5 cm and contains 
various services, including cooling, power cables, and silicon detector readout. The crack is 
not projective to the interaction point. In addition, about 1% of the ECAL crystals are either not 
operational or have a high level of electronic noise [17|, and they are masked in reconstruction. 
The rj-(p distribution of these crystals for the barrel, and y-z distribution for the endcaps, is 
shown in Fig. [3] 

In this Section, we illustrate the effect of these features on the Bj distribution and test the 
reliability of the simulation for events with jets pointing towards masked ECAL channels or 
cracks. 

Figure |4] shows the B j distributions from simulated samples of events containing at least 2 
jets, with the leading jet satisfying pj > 50 GeV and the second jet satisfying pj > 25 GeV, 
for Calo^i/ TCKt/ and PF£t- For events with 100 < Bi < 200 GeV, the contribution from 
QCD multijet production is 24-42%, depending on the B j reconstruction algorithm; the rest is 
from W/Z/tt production. In order to illustrate the effect of the cracks, distributions are also 
shown for those subsets of these samples that have at least one jet that is aligned with the Bj 
to within A</>(£x/jet) < 0.2 and that is pointing towards masked ECAL channels, the barrel- 
endcap boundary (1.3 < \rj\ < 1.7), or the endcap-forward boundary (2.8 < \rj\ < 3.2). The 
masked ECAL channels considered here are those that are part of a group of 5 x 5 or 5 x 1 
masked channels that are adjacent in rj-cp, as they have larger impact on the Bj distributions 
than isolated masked channels. A jet is considered to be pointing to one of the masked ECAL 
channels when its jet centroid is within AR < 0.2 of a masked ECAL channel, where AR = 
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Figure 2: (left) Probability of finding a beam-halo tagged event in muon-triggered events. Re- 
sults are shown as a function of the beam intensity, (right) PF B j distribution for all the events 
from muon and Calo B j triggers that were analyzed, and for the subset of these events that 
were identified as beam halo. 
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Figure 3: Distribution of masked ECAL channels in (top) barrel: t]-(p view of 170x360 individual 
crystals and (bottom) endcaps: y-z view of 2x7400 individual crystals. 

y(A77) 2 + (A</>) 2 . We can see effects on the B j distribution from masked ECAL channels, while 
the calorimeter boundaries do not appear to have an enhanced contribution to the events with 
large Bj. 

Figure [5] shows the fraction of dijet events with at least one jet aligned with the B j and also 
pointing towards the masked ECAL channels, the barrel-endcap boundary, or endcap-forward 
boundary for data and for simulation. Figure |5peft) shows that the masked ECAL channels 
enhance the rate of events with large B j in both data and Monte Carlo simulation. Approx- 
imately 20% of the events with Bt > 80 GeV have contributions to the measured Bt from 
mismeasurements due to masked ECAL channels. Results from simulations indicate that the 
fraction of events with large B t due to mismeasurements (excluding the predicted contribu- 
tions from sources of genuine Bt, such as W/Z/tt) is 30%. As shown in Figs. |5|middle) and 
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area, for events containing at least 2 jets with p^ 1 ^ > 50 (25) GeV (black solid), and for the 
subsets of these events with a jet aligned with Uj within A</>(£x/jet) < 0.2 and pointing to- 
wards a masked ECAL cell (red dotted), the barrel-endcap boundary (green dot-dashed), and 
the endcap-forward boundary (blue dashed) in simulation. 
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Figure 5: Fraction of dijet events in data (points) and simulation (red band) with a jet aligned to 
Kt within A^Kivjet) < 0.2 and pointing towards (left) a masked ECAL channel, (middle) the 
barrel-endcap boundary, and (right) the endcap-forward boundary, in data and in simulation. 

fright), the fraction of events which contain a jet that is both aligned with the $j and pointing 
towards a calorimeter boundary does not have a strong dependence on $ j, indicating that the 
calorimeter boundaries are not major contributors to events that have large apparent B j due 
to mismeasurements. Unlike the masked ECAL channels, the cracks are not projective to the 
interaction point, and therefore energies of particles traversing these cracks are still measured, 
albeit with degraded resolution. 

While the impact of the cracks is small, analyses sensitive to events with large # j need to take 
the ECAL masked channels into account. About 70% of the ECAL channels that are masked 
during offline reconstruction have a useful measurement of their energy from the separate 
readout of the LI trigger. Although the trigger readout saturates, it can be used to recover 
energies smaller than this and to identify events that had more than this amount of energy in 
a masked channel. This saturation energy has been increased from 64 GeV to 128 GeV in 2011. 
Analysts can veto events with a jet pointing towards an ECAL masked channel that does not 
have trigger information or that has trigger-readout energy at the saturation threshold. 
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6 Missing transverse energy scale and resolution 

In this section, we study the performance of B j using events where an identified Z boson or 
isolated 7 is present. Events containing vector bosons may be produced in hard parton-parton 
collisions such as qg— > q7, qq — > Z, qg— >qZ, and qq — > gZ. While there is no genuine It in 
these events, we can induce it by removing the vector boson. By comparing the momenta of 
the well-measured and well-understood vector boson to the B j induced this way, we probe the 
detector response to the global hadronic system and measure the scale and resolution of Bj- 
While the lowest order underlying processes may be simple, many physics and experimental 
issues contribute to the measured, induced Bj in these events. Effects due to jet energy scale 
corrections and fluctuating jet composition directly impact the measurement of the hadronic 
products of the hard collision. Underlying event activity, pile-up, detector noise, and finite 
detector acceptance contribute as well. 

The following notation is used: the vector boson momentum in the transverse plane is qj, 
and the hadronic recoil, defined as the vector sum of the transverse momenta of all particles 
except the vector boson (or its decay products, in the case of Z candidates), is Uj. Momentum 
conservation in the transverse plane requires qj + Uj = 0. The recoil is the negative of the 
induced Bt- 

The presence of a well-measured Z or 7 provides both a momentum scale, qj = \qj\, and a 
unique event axis, qj. The hadronic recoil can be projected onto this axis, yielding two signed 
components, parallel (uy) and perpendicular (u±) to the event axis. Since Un = uj ■ cjj, and the 
observed hadronic system is usually in the opposite hemisphere from the boson, Uy is typically 
negative. 

The mean value of the scalar quantity (u^)/qj is the scale factor correction required for Bj 
measurements in the classes of events considered here, and is closely related to jet energy scale 
corrections and jet parton flavour. We refer to (u^)/qj as the "response" and denote distribu- 
tions of this quantity versus qj as "response curves". Deviations of the response curve from 
unity probe the Bj response as a function of qj. 

Resolution is assessed by measuring the RMS spread of My and about their mean values, 
after correcting for the response, and is denoted RMS(wy ) and RMS(u±). As with the response, 
we examine the resolutions as functions of qj. 

6.1 Direct photon sample 

Candidate photon events are selected by requiring each event to contain exactly one recon- 
structed photon in the barrel portion of the ECAL (\t]\ < 1.479), with qj > 20 GeV, and which 
passes the identification and isolation selection described in Section [3] The total number of 
events passing all requirements is 15 7567, of which 67621 have only one reconstructed pri- 
mary vertex. The prescale factors for the HLT triggers used to collect this sample varied over 
the course of the 2010 LHC running period. As a result, this sample is dominated by events 
recorded during the earlier period of the data taking, when the fraction of crossings containing 
pile-up interactions was smaller. 

Figure [6] shows the photon qj spectrum for data and for simulation. About half of the observed 
rate arises from QCD dijet production where one jet passes all photon identification require- 
ments. Such jets are typically highly enriched in 7Z° — > 77 and contain little hadronic activity. 
The detector response to these jets is similar to that of single photons, and studies indicate 
that response curves extracted from these QCD background events match the response of true 
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6 Missing transverse energy scale and resolution 
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Figure 6: Distribution of qj for events selected as photon+jet candidates. Predicted rates from 
simulation for signal and backgrounds are also shown. QCD refers to multijet events. 

photon-jet events to within a percent. We therefore make no further attempt to filter them out. 
The detector response to the jet depends on the type of parton from which it originated. The 
leading jet in photon events is predicted to predominantly be a quark jet. A prediction for the 
difference in response for the CMS detector between quark and gluon jets can be found in [20 1. 
The difference is largest for Calo Jets (~ 20% for jets with pj of 20 GeV), and decreases with pj. 
The primary reason that the response is lower for gluon jets is that their particles tend to have 
lower pjs, and the calorimeter response is lower at low pj. For PF jets and PF Ktv which use 
tracker instead of calorimeter momenta for most charged hadrons, the difference in response is 
reduced, and varies from about 5% at 20 GeV to a percent a high pj. 

6.2 Z samples 

For the Z — > e + e~ selection (electron channel), we require two well-identified and isolated 
electrons, with pj > 20 GeV, within the fiducial region of the ECAL. The invariant mass (M«) 
of the electron pair is required to be in the range 70 < M« < 120 GeV. 

For the Z — > selection (muon channel), we require two isolated muons with opposite 

electric charges, that have pj > 20 GeV, and are within the < 2.1 region. The invariant mass 
Mm of the muon pair is required to be at least 60 GeV, and no more than 120 GeV. 

We obtain a total of 12 635 (12383) Z — > e + e (Z — > n + }T) candidates. The relative contribu- 
tions of signal and background are estimated from simulation. By normalizing the invariant 
mass distribution from simulation of signal and background so that it has the same number of 
events that is observed in the data, a total background of around 143 (35) events is estimated 
to have a contribution of 97 (2) events from QCD, 28 (9) events from electroweak, and 18 (24) 
from final states containing top quarks. 

Figure [7] shows the Meg distribution for the electron and muon samples. Figure [8] shows their 
cjj spectrum. Except at very low qj, the leading jet in Z events, as with the 7 events, should 
usually be a quark jet. 



6.3 Scale and resolution for events with one primary vertex 
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Figure 7: (left) Invariant mass distribution of the two leading electrons and (right) invariant 
mass distribution of the two leading muons, for the Z boson candidates, along with the pre- 
dicted distribution from simulation. QCD refers to multijet events. 
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Figure 8: The qj distribution for Z boson candidates in (left) the electron channel and (right) the 
muon channel, along with the prediction from simulation. Systematic uncertainties are shown 
as grey bands. 

6.3 Scale and resolution for events with one primary vertex 

To study the Bt scale and resolution, we decompose the recoil with respect to the boson (7 
or Z) direction in the transverse plane. We restrict ourselves to events containing one recon- 
structed primary vertex. Z yields are thus reduced to 2611 Z — > e + e candidates and 2438 
Z — > candidates. The effect of pile-up on the scale and the resolution is studied in Sec- 

tion 6.5.1 Distributions of the components of the recoil calculated from PF Bj that are parallel 
and perpendicular to the boson axis, Un and u±, are shown in Fig. [9| for direct photon candi- 
dates, Z — > e + e~ candidates, and Z — > candidates. As expectea, the parallel component 
is mainly negative, consistent with the back-to-back nature of the events, while the perpendic- 
ular component is symmetric. 



The distributions are corrected for the residual contamination (5 ± 1%) from events with more 
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Figure 9: U\\ distributions for PF^x for (top, left) 7, (top, middle) Z — > e + e~ and (top, right) 
Z — > events; distributions for PFBj for (bottom, left) 7, (bottom, middle) Z — > e + e~, 

and (bottom, right) Z — > events. Distributions are for events containing one PV, with 

a correction for residual pile-up contamination. Statistical uncertainties and systematic uncer- 
tainties on the removal of events with more-than-one PV and on the boson qj spectrum are 
shown as grey bands on the prediction from simulation. QCD refers to multijet production. 

than one interaction. The number of events with more than one scattering that are recon- 
structed as a single PV is estimated by convoluting the efficiency for reconstructing two vertices 
as a function of the vertex separation with the z distribution of vertices. The distributions are 
corrected for this contamination by subtracting multi-vertex-event shapes obtained from data, 
rescaled to the estimated contamination, from the distribution from events with one PV. The 
systematic uncertainty on the residual contamination is obtained by varying the normalization 
within its uncertainties. 

Events generated with pythia are reweighted so that the qj spectrum matches that predicted 
by the ResBos Monte Carlo program |2TI , in order to take advantage of its resummed calcula- 
tion of the boson qj spectrum. The systematic uncertainties due to our imperfect knowledge of 
the true qj distributions for Z bosons are estimated from the difference between the qj distri- 
butions predicted by pythia and ResBos. We set the systematic uncertainty, bin-by-bin in qj, 
equal to this difference. 

In addition, there is a systematic uncertainty on the prediction from the simulation due to 
the size of the simulation samples. In Fig. |9j the dominant uncertainty on u\\ for the part of 
the distribution from the electron channel with u\\ < — lOOGeV is from this source and from 
uncertainties on the removal of the multi-PV contamination. 



Figure 10 shows the response curves, |(mi|)|/^t versus qj, extracted from data, for the three 
£t reconstruction algorithms, Calo£x, TCKt/ and PFKt for 7, Z — > e + e , and Z — > 
samples. The agreement in response between the different samples is good. The agreement 
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Figure 10: Response curves for events with one primary vertex, for (left) Calo^x, (middle) 
TC£j/ and (right) PF£x- Results are shown for photon events (full blue circles), Z — > e + e~ 
events (open red circles) and Z — > \i^}T events (open green squares). The upper frame of each 
figure shows the response in data; the lower frame shows the ratio of data to simulation. The 
vertical axis labels at the far left apply to all three subfigures. 

between data and simulation is good, and the results indicate that the three reconstruction 
algorithms are distinct in their capabilities, performing differently in the recovery of hadronic 
activity in the detector. The response for Calo£j is slightly larger than one because the jet 
energy scale used in the type-I corrections was determined from a sample with a mixture of 
quark and gluon jets, while for these samples the leading jet is primarily a quark jet. The TC B j 
response is lower because it has neither type-I nor type-II corrections. The PF K t response is 
lower than the Calo B j response at low values of qj because Calo £ x has type-II corrections 
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Figure 11: Resolution curves for components of hadronic recoil measured in direct photon can- 
didate events with one primary vertex, (left) parallel to boson; (right) perpendicular to boson. 
Data and simulation are indicated by points and histograms, respectively. Black circles (up- 
per): Calo Ex; Pink squares (middle): TCKx; Blue triangles (bottom): PF^x- Shaded regions 
indicate statistical uncertainties on the simulation. 
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6 Missing transverse energy scale and resolution 
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Figure 12: Resolution curves for components of hadronic recoil calculated using PF Ex mea- 
sured in events with one primary vertex. Left: parallel to boson; right: perpendicular to boson. 
Z — > , Z — > e + e and 7 events are indicated by open green squares, open red circles and 
full blue circles, respectively. The lower frame indicates the ratio of data to simulation. The 
vertical axis labels at the far left apply to both subfigures. 

while PF £x has only type-I corrections. 

Figure [TT]shows the resolution curves from photon candidate events for un and for data and 



simulation, for all three reconstruction algorithms. Figure 12 shows the resolution as measured 
in 7, Z — > e + e~, and Z -4 events for PF The measured resolution must be corrected 

for the scale to avoid a misleading result; e.g., the apparent resolution on u± is proportional to 
the scale and therefore an algorithm with a scale that is smaller than unity could appear to have 
a better resolution than one with a scale of unity without such a correction. Since only Calo M t 
has been corrected fully for the detector response with both type-I and type-II corrections, the 
resolution measurements are rescaled, bin by bin, using the corresponding response curves of 



Fig. 10 The data confirm the prediction from simulation that tracking information significantly 
enhances the £ x resolution. The resolutions as measured in the different samples are in good 
agreement, but are ~ 10% worse than expected from the simulation. A similar difference in 
resolution for jets for the 2010 run is documented in [14|. The small discrepancies between data 
and simulation shown in Fig. [9] are due to this difference. 



6.4 Resolution in multijet events 

The Kt resolution can also be evaluated in events with a purely hadronic final state, where 
the observed B j arises solely from resolution effects. Because the $ j resolution has a strong 
dependence on the associated £ Ej, it is presented as a function of £ Ej. We characterize the 
Bt resolution using the a of a Gaussian fit to the distribution of the x and y components of 
fx (Ex,y)- In order to make a meaningful comparison, we calibrate the measured Kx for the 



different algorithms to the same scale using the response from Fig. 10 These corrections would 
not be needed if all types of Uj had both type-I and type-II corrections. 



6.5 Effect of multiple interactions 
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Figure 13: Calibrated t x ,y resolution versus calibrated PF £ Ej for Calo £x, TC £x, and PF £x 
in data and in simulation. 



For £ Ej, we use the PF £ £j as measured by the particle-flow algorithm for all types of £ j, as it 
gives the best estimate of the true £ Ej, and hence is an accurate evaluation of the event activity. 
We use PF Ej for all algorithms to ensure their measure is the same. We calibrate PF Ej to 
the particle-level £ Ej, on average, using the predicted average mean value as a function of the 
particle-level Ej from a simulation of events from the pythia 8 event generator 



Figure 13 shows the calibrated Zx,y Gaussian core resolution versus the calibrated PF TJEx for 
different £ t reconstruction algorithms in events containing at least two jets with pj > 25 GeV. 
Both TC £x and PF £x show improvements in the £x resolution compared to the Calo Bj, and 
the PF % x yields the smallest £ j resolution. 



Figure 14 shows the PF £ t distributions for different intervals of Calo £j and for jet mul- 
tiplicities varying from two to four, normalized to the same area. The jets are required to be 

in- 



above a pj-threshold Q f 20 GeV. The good agreement of the normalized shapes in Fig. 14 
dicates that PF £ x-performance in events without genuine £ x is driven by the total amount of 
calorimetric activity (parametrized by Calo £x) and no residual non-linear contribution from 
jets to PF £x is visible. Similar behaviour is also observed for Calo £x and TC £x- 

6.5 Effect of multiple interactions 

Pile-up, namely multiple proton collisions within the same bunch crossing, occurs because of 
high LHC bunch currents and can play an important role in £x performance. 

Because there is no true £ x in minimum bias events and because the average value for a com- 
ponent of £x in these events is zero (e.g., the x or y component), pile-up should have only a 
small effect on the scale of the component of the measured £x projected along the true £x di- 
rection. Pile-up, however, will have a considerable effect on the resolution of the parallel and 
perpendicular components. 



We investigate the effect of pile-up using multijet samples, 7, and Z data. 
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6 Missing transverse energy scale and resolution 
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Figure 14: PF Bt distributions in 2-, 3- and 4-jet events, in selected Calo £ Ej bins. 
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Figure 15: The ratio of the response for the component of the induced along the boson 
direction, measured in 7 events for events containing 1 PV and at least 2 PVs for (left) Calo Kt, 
(middle) TC Bt, and (right) PF % j. Also given is the best fit value for the average ratio, which 
corresponds to the solid, red line. 

6.5.1 Studies of pile-up effects using photon and Z events 

In this section, we use samples containing a vector boson to measure the effect of pile-up on 



the scale and resolution of a component of $j. Figure 15 shows the ratio of the response as 
measured in 7 events for events containing 1 PV and at least 2 PVs. The ratio of the responses 
of the component of the measured Bt along the boson direction is close to one, as expected. 
It is slightly larger at low qj when pile-up is present. This is expected, as pile-up can reduce 
energy lost due to zero suppression in the readout of the calorimeter if energy from a pile-up 
interaction and from the hard scattering are both in the same readout channel. If the sum is 
larger than the zero suppression thresholds, more of the energy from the hard scattering is 
recorded. 



Figure 16 shows the resolution versus the of the 7 for the components of the hadronic recoil 



parallel and perpendicular to the boson direction for 1, 2, and 3 reconstructed PVs. Also shown 



is the prediction from simulated 7 events without pile-up. Figure 17 shows the resolution 
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Figure 16: Resolution versus photon qj for the parallel component (top) and perpendicular 
component (bottom) for (left to right) Calo B j, TC B j, and PF B x, for events with 1 (circles), 2 
(squares), and 3 (triangles) reconstructed primary vertices. 

versus the qj of the Z for the parallel and perpendicular components of the hadronic recoil. 



The parameterization of Bj resolution used in Figs. 16 and 17 is given by : 



^otal = (« V 7 ^ + b ) 2 + Koise/ ES (<?T)) 2 + (N - 1) ( WesM) 2 (1) 

where a and b characterize the hard process, cr noise is the intrinsic noise resolution, N is the 
number of reconstructed vertices in the event, cr PT j is the intrinsic pile-up resolution, and /es( ( ?t) 
is the energy scale correction applied on each event. At low qj, the resolution is dominated by 
contributions from the underlying event and detector noise (c no j se ). Since these contributions 
can not be distinguished from those due to the particles from the recoil, and since the recoil 
measurement needs to be corrected for the detector response, these contributions are magnified 
and have a larger contribution at low boson qj when energy scale corrections are applied. As 
expected, the resolution is degraded with increasing pile-up interactions. Results from the Z 



and 7 channels are in agreement and are similar to the values obtained in Section 6.5.2 from jet 
data. 

6.5.2 Studies of pile-up effects in jet data 

In this Section, we study the behaviour of the PF B t distributions in samples containing high pj 
jets when pile-up is present. The data are selected using a prescaled Hj trigger with a threshold 
of 100 GeV, where Hj is defined as the scalar sum of the transverse momenta of PF jets (pj > 20 
GeV, \rj | <3). Additionally, in the offline analysis, each event is required to have Hj (calculated 



using PF jets) > 200 GeV to avoid bias from the trigger. Figure 18 shows that the widening of 



the PF £ x distribution with increasing number of vertices can be modeled by convoluting the x- 
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Figure 17: Resolution versus the qj of the Z for the parallel component (top) and perpendicular 
component (bottom) for (left to right) Calo B j, TC B j, and PF B j, for events with 1 (circles), 2 
(squares), 3 (triangles), and 4 (diamonds) reconstructed primary vertices. 

and y-components of the one-vertex Bt shape with a Gaussian (G) whose mean is (n — 1) • Ay. x 
and with standard deviation \/n — 1 ■ Aa x : 



lj,n = V (2*1 ® G[(n - 1) • A^ X/ Vn - 1 • A^]) 2 + (l yl ® G[(n - 1) • A^ y , Vn - 1 • A^]) 2 

(2) 

where B x ,y are the x and y components of £j- Here we assume that each additional vertex 
contributes with a constant Aa x (Act,/) to the B j resolution such that the resolution with n pile- 
up interactions is related to that with one primary vertex via: a xn = cr^i + {n — l)Ac 2 . In 
addition we also allow for a linear shift of B x and By by A}i x (A}iy) such that ]i xn = ]i x \ + {n — 
l)Aji x . A fit of Eq. (|2| to data results in Aa x = Aa y = 3.7 GeV, consistent with the results from 



Section 6.5.1 This fit is performed simultaneously on the Bj distributions of events containing 
two to seven vertices. The shifts of the x and y Bj components are estimated to be Aji x = 0.5 
GeV and A]i XJ = —0.3 GeV respectively, which are small compared to Ac and are consistent 
with the expected shift seen in simulation due to non-functional channels. 
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Figure 18: PF£x distributions in pile-up events. The figures show a comparison between the 
one- and n-vertex shapes (n = 2.. .7) and the results of a simultaneous fit of Eq. (|2j to the n-vertex 
shapes. The one-vertex distribution is normalized to the n-vertex distribution for each plot. 
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7 Studies of physics processes containing genuine £ T 

In this section, we examine distributions relevant to £x in events containing neutrinos. Events 
containing W bosons and b quarks are studied. 

7.1 W events 

The performance of £x is studied in events that contain large, genuine Uj- W — > £v events, 
where £ is a muon or electron. For most W events, the magnitude of £ x is approximately equal 
to the px of the charged lepton, but its resolution is dominated by the hadronic recoil. When 
the W cjj is small compared to the W mass, the £ t is approximately 



£t ~ PtCO -0.5ui 



where is the component of the recoil parallel to the lepton transverse direction. 

In the W — > ]iv decay channel, events are required to have been collected by a single-muon 
high-level trigger. In addition, candidates are selected by requiring a muon with \rj\ < 2.1 that 
has px > 25 GeV. Events with a second muon with pj > 25 GeV are rejected to suppress Z 
and tt contamination. W — > ev decays are identified using similar selection criteria. A single- 
electron high-level trigger requirement with a px threshold of 15 GeV is applied. Events are 
also required to contain an electron with pj > 25 GeV. Events with a second electron with 
pj > 20 GeV are rejected, and rejection against 7 conversions is applied. A total of 24628 
(29 200) W — > jiv (W — > ev) events with only one primary vertex are selected. 

The main sources of background to the W — > £v signal are jet events with one jet falsely iden- 
tified as a high-px muon or electron and Z — > ££ events with one lepton escaping detection. 
The jet events usually have low values of £x- The apparent £x in these events (that have no 
genuine £ x) is amplified by scale corrections to the £ x (type-I and type-II) because, since they 
are indistinguishable, artificial £ x receives the same scale correction factor as genuine £ x, and 
genuine £ j tends to be underestimated without corrections. Other backgrounds include W and 
Z bosons decaying into T, followed by T — > Ivv, and tt events, with one top quark decaying 
semi-leptonically. The relative normalization of the different electroweak (EWK) signal simula- 
tion event ensembles, and simulations of those backgrounds that contain an electroweak boson 
(W — > £v, Z — > ££, tt), are set by the ratios of their theoretical cross sections computed at next- 
to-leading order [231. The normalization of the composite EWK and the QCD contributions are 
established through a one-parameter binned fit to the £ x distribution from data. 



Figure 19 shows the PF £x distribution for the W — > ev and W — > }iv candidate samples, along 
with the expectation from simulation. As for the analyses using Z events, the background 
distributions include a grey-shaded band indicating the estimated uncertainty due to the size of 
the simulation samples, modeling of the W qj spectrum, and the pile-up correction procedure. 
In most cases this uncertainty is too small to be visible. Data and simulation agree well, and 
the W shows up prominently as expected. 



Figure 19 also shows the uj distribution. To suppress QCD background, for the uj studies 
only, we further require that the W candidate pass a requirement on the transverse mass, de- 
fined as Mx = ^/ pj(£) ■ £x ■ (1 — cosA<^), where A<p is the opening angle in the transverse 
plane between the lepton candidate and the £x- We require Mj > 50 GeV, and a minimum 
£x threshold, £x > 25 GeV. The £x resolution has substantial contributions from the mis- 
measurement of the many particles in the underlying event. These contributions can be more 



7.2 Heavy flavour production and E T 



23 



clearly seen in the iij distribution, since they are not obscured by the contributions from the 
charged lepton. Again we see reasonable agreement between data and simulation. 



.3000 




50 100 
PF g T [GeV] 



3000 



CD 
> 
CD 



: 2000- 



1000 




40 60 
PF [GeV] 



> 

CD 

a 



CD 
> 
CD 



CD 

E 

=3 



x10 



36 pb" 1 at\fs = 7 TeV cms 



• Data 

■ Background 




60 80 100 
PF £ T [GeV] 



> 4 

CD 

o 

CM 



CD 
> 
CD 



x10 



CD 



□ 



36 pb" 1 at\/s = 7 TeV cms 

• Data 
| | W -> n v 
■ Background 




40 60 
PF u T [GeV] 



Figure 19: The PF t t (left) and |mt| (right) distribution in W — > ev (top) and W — > ]iv (bottom) 
candidate events. Both data (points) and simulation (solid lines) are shown. The plots on the 
right include selection requirements on £ j and Mj, while those on the left do not. 



7.2 Heavy flavour production and £ T 

The Bt distributions from jet samples containing b quarks can differ from those of inclusive 
jet samples because the B hadrons have unique fragmentation properties, and sometimes their 
final states contain neutrinos. Neutrinos from b jets are one of the main sources of severe 
underestimations of jet energies (the other main source is the ECAL masked channels discussed 
in Section [53)l . In this section, we study the induced £j in an inclusive b-tagged jet sample. 



We compare the Bt distributions in dijet events with and without a secondary vertex, i.e., 
events where the leading jet has a positive SSV (displaced vertex) tag or a SMbyPt (lepton) 
tag 00 (Section [3] includes a description of these algorithms). Since these taggers rely on 
tracking information, we require that the leading jet has \rj\ < 2.1. Also, we require the lead- 
ing two jets to have pj > 40 GeV. Below this value, the b-quark tagging purity is significantly 
reduced. A prescaled jet trigger with a p x threshold of 15 GeV was used; the resulting sample 
corresponds to an integrated luminosity of 0.025 pb -1 . Figure 
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shows the fraction of events 
from this sample with a b-tagged jet as a function of £ j for the two tagging algorithms. The 
larger increase in the fraction of b-tagged events at large U j for SMbyPt than for SSV is due, in 
part, to the higher probability of neutrinos in leptonically tagged events. 
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Figure 20: Fraction of events in a jet sample that contain a jet tagged as a b jet by the SSV- and 
SMbyPt-taggers for (left) Calo t t and (right) PF I T . 

8 M T significance 

A spurious nonzero £ j in an event can have contributions from many sources, including mea- 
surement resolution, reconstruction inefficiencies, instrumental defects, and improper pattern 
recognition. Events in which the reconstructed £ t is consistent with contributions solely from 
particle-measurement resolutions and efficiencies can be identified by evaluating the £t sig- 
nificance, S. The significance offers an event-by-event assessment of the likelihood that the 
observed £t is consistent with zero given the reconstructed content of the event and known 
measurement resolutions. A similar variable used by the CDF collaboration is described in 



8.1 Definition 

The significance requires evaluation of the uncertainty in the total measured transverse energy, 
which is given by 



eT 31 = E h 

iex 



-t- 



(3) 



where Ej. = (E Xi , E yi ) is the measured transverse momentum of the z' th reconstructed object. X 
is the set of reconstructed objects, such as calorimeter towers (for Calo Ex) or PF particles (for 
PF Et)/ used to calculate Et- In the derivation of the significance, there are three relevant quan- 
tities for each object in the sum. The first of these is e^, the true transverse momentum of the 
object. The significance provides a measure of whether an event is consistent with the null hy- 
pothesis of zero genuine total transverse momentum. Under this hypothesis, Ya ex^Ti = 0- The 
second is Ej., the measured transverse momentum of the object, which is distributed accord- 
ing to Pi{Ej i \ej.J, the probability density function (pdf) for observing the measured transverse 
momentum given the true transverse momentum of the object. The third is e ; = Ej. — ej r For 
convenience, we define an equivalent pdf in terms of this difference: (e, \ej.) = Pj (tfj + ej. \ej { ). 
Given the null hypothesis, TJEii = E^v so that the z' th reconstructed object contributes £, to the 
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measured total transverse momentum. 

We first introduce the likelihood that we would observe a total transverse momentum e under 
our null hypothesis. For the two object case, the likelihood function is given by 

£(e) = |p 1 (E T Je T J)P 2 (E T J%))^(e-(E Tl +E T2 ))rfE T2 dE T2 
= J p 1 (e 1 \e ll )p 2 (s2\ej 2 )5(e- (e 1 + e Tl +£2 + eT 2 ))d£ 1 de 2 
= J pi(e 1 \e ll )p2(e 2 \eT 2 )S(£- (e x + £ 2 ))de 1 de 2 , (4) 

since = Ya^i = <h\ + ^t 2 - F° r an arbitrary number of input objects, the full likelihood 
function can be generated by a recursive application of Eq. Q. The significance is defined as 
the log-likelihood ratio 

which compares the likelihood of measuring the total observed Ej — D Ejj =1] £j to the 
likelihood of the null hypothesis, Ex* = 0. 

This formulation is completely general and accommodates any probability distribution func- 
tion. In practice, however, we often employ Gaussian uncertainties for measured quantities, 
for which the integrals of Eq. (|4j can be done analytically. The Gaussian probability density 
function is given by 



Pi{ti\e Ti ) ~ exp (-^(t i ) T V i 1 (e i ^J , 



where V, is the 2x2 covariance matrix associated with the z th measurement. The integration 
of Eq. Q yields 

with V = Vi + V2. When many measurements contribute, the expression generalizes to 



^)~exp(-i(e) T ^£v^ (e)V 



(6) 



The covariance matrix U; for each reconstructed object in the Ej sum is initially specified 
in a natural coordinate system having one axis aligned with the measured Ejj vector, £7. = 
(E Ti cos (pi, E T . sincpi): 
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(We adopt the simplifying assumption that Ej and <p measurements are uncorrected.) This 
matrix is rotated into the standard CMS x — y reference frame to give the error matrix 

V i = R(<p i )U i R- 1 (<p i ), (8) 

where R (</>,) is the rotation matrix. The 5 -matrix summation is then performed in this common 
reference frame. Combining Eqs. ([5]), (|6|, and (|8j yields 





U/R- 1 (</>;) L £ t, • (9) 



Equation (j^Jl makes explicit the dependence of S and £x on the set of objects X over which 
the vectors and matrices are summed. In general S is small when the B j can be attributed to 
measurement resolution, and large otherwise. 

In the Gaussian case, S is simply a x 2 with two degrees of freedom. If we rotate into a coor- 
dinate system with the x axis parallel to the B j axis, instead of the CMS horizontal axis, then 
Eq. (9| is simplified to S = E\l (c| T (l — p 2 )), where is the variance of the magnitude of Bj, 
and p is the correlation coefficient between the variances parallel to and perpendicular to the 
measured Bj- This form emphasizes the essential meaning of S, but obscures the important 
feature that, through its denominator, S embodies the full topological information in the event. 
Essential features such as the angles between the measured B t arid the reconstructed objects 
in the event are embedded in the definition of the denominator. This form also makes appar- 
ent the relationship between the true significance (in the Gaussian limit) and the more naive 
measure E = BtI VX)Et- 

The specialization to a Gaussian probability density function is less restrictive than it may ap- 
pear, as any probability density function expressible as a linear combination of Gaussians is 
accommodated by the formalism presented here. 



To apply Eq. (j9| to PF Bj significance, we note that the Gaussian pdf only accommodates mea- 
surement resolution. Using only reconstructed PF particles to determine the covariance matrix 
would neglect fluctuations in the measured PF particle content itself. These fluctuations arise 
from finite detection and reconstruction efficiencies, and provide a non-negligible contribution 
to the PF £ x resolution. These fluctuations, however, also affect the PF jet resolutions. We can 
therefore substitute the PF jet resolutions for the combined measurement resolutions of the PF 
particles that have been clustered into jets. Hence the sum of covariance matrices in Eq. Q in- 
cludes contributions from PF jets, PF particles that were not considered during jet finding (e.g. 
isolated leptons), and PF particles that are not clustered into any jet. This approach inherently 
takes into account the contributions both from measurement resolution and from fluctuations 
in the reconstructed particle content. 

The covariance matrices U, of Eq. |7j| are obtained from our knowledge of the response of each 
type of PF particle or jet as a function of pj and tj. The charged hadron and muon resolutions 
are obtained on a particle-by-particle basis from the error matrix from the final track fit, and the 
resolutions for electrons are those obtained from studies of data samples of known resonances 
such as neutral pions, the Z, etc. The jet and photon resolutions are from simulation. No input 
resolutions were tuned based on the behaviour of the significance distribution itself. 
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8.2 Performance of <S PF in dijet events 

Because 5pf is ^ 2 distributed, it should exhibit a flat probability of x 2 , ^P{x 2 )r f° r two degrees 
of freedom in an event sample that nominally has no genuine £j. (That is, 1 — V(x 2 ) is the 
standard cumulative distribution function of the x 2 statistic for two degrees of freedom.) Dijet 
samples from pp collisions are dominated by such events. 

We select dijet events by requiring at least two jets satisfying \t]\ < 2.3 and pj > p™ in , with 
thresholds p™ m of 30 or 60 GeV. One of the jets above threshold must have been responsible 
for the event passing an HLT single-jet trigger. We use data collected with a 15 GeV trigger 
threshold for our 30 GeV dijet sample, and a 30 GeV trigger threshold for the 60 GeV dijet 
sample (Because of different prescale factors applied to the two trigger streams, the 60 GeV 
dijet sample is not a direct subset of the 30 GeV sample.) 

We compare the 5pp distributions, as well as their corresponding V{x 2 ) distributions, in data 



and simulation in Fig. 21 for both values of the p™ m threshold. The significance distribution 
very closely follows a pure exponential, and the V{x 2 ) distribution is populated quite uni- 
formly between zero and unity in both data and simulation. There is a small peak at zero in 
V(x 2 )', simulation (Fig.|2l]l indicates that about half of this peak results from genuine £x in the 
event sample. This Bt arises from a combination of sources such as the semileptonic decays of 
heavy quarks and the finite t] acceptance of the detector. The data and simulation distributions 
match well in the 30 GeV threshold sample. MC studies show that the remainder of the excess 
of low probability events after accounting for genuine £j typically have at least one high- p j 
jet whose response is in the non-Gaussian tail of the response function. 

To probe the stability of the 5pp behaviour, we have studied dijet samples with different p™ m 
thresholds, which changes the relative contributions of different detector regions in the covari- 
ance matrix calculations. We find that, overall, the 5pp distributions for the bulk of the data 
continues to exhibit near-ideal behaviour independent of threshold. As the 60 GeV sample 
shown here demonstrates, though, the higher threshold data does begins to develop a larger 
tail in the significance, and a correspondingly larger peak at zero in V(x 2 )' than we find in the 
simulation. The discrepancy between data and MC is below the 0.2% level. Visual examina- 
tion of the events with low probability reveal that the discrepancy arises from a combination 
of events with a residual anomalous energy contamination and other events with a high-px 
jet with activity straddling the endcap (HE) and forward (HF) calorimetry for which the non- 
Gaussian tails are not yet perfectly modeled. 

For the 5pp distributions shown here, the transition point for use of resolutions based on PF 
jets rather than resolutions from unclustered PF particles in the 5 P f calculation (Eq. Q) occurs 
at a jet pj of 3 GeV. The 5pp distributions are insensitive to the variation of this transition point 
between jets and individual particles over the range of 1 to approximately 6 GeV. By 10 GeV, 
a slope in the V(x 2 ) distribution has clearly appeared, indicating that we no longer account 
sufficiently for contributions to the £ t resolution from fluctuations in the reconstructed particle 
content. 

A powerful feature of the £t significance is that its distribution is insensitive to pile-up (for 
events with no genuine Bt)- As long as the correct resolutions are input, the significance should 
still have a pure exponential behaviour with a uniformly distributed V(x 2 )- In Fig. 
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no 



restrictions were made on the number of interaction vertices in the data, while the simulation 



has no pile-up. In Fig. 22 we compare the shapes of the single-vertex and multiple-vertex 
significance and V(x 2 ) distributions in data. The shapes are very similar, as expected. The main 
difference arises in the low probability region, where the multiple interaction data exhibits 
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Figure 21: The £j significance tSpp distributions (left) and the corresponding probability of 
X 2 , V{x 2 ), distributions (right) for dijet event samples in data (points) and simulation (solid 
histograms) with 30 GeV (top) and 60 GeV (bottom) jet pj thresholds. The dashed histograms 
show the simulation distributions with true Bt contributions, from physics and finite accep- 
tance effects, subtracted event-by-event. The dotted line overlaid on the 5pf distributions 
shows a reference pure exponential function. Each inset expands the small V(x 2 ) region. 
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behaviour closer to the ideal - an example of the central limit theorem. With the additional 
contributions to the Bj resolution, the roles of the non-Gaussian response tails and genuine 
Bj are diminished. The overall insensitivity can be useful, for example, when extrapolating 
backgrounds dominated by samples with nominally zero genuine B j. 




Figure 22: The Bj significance (left) and V(x 2 ) (right) distributions for events with a single 
interaction vertex (histogram) and multiple interaction vertices (points) in the 60 GeV threshold 
dijet sample. The inset expands the small V(x 2 ) region. 

8.3 Application to W — > ev events 

As a case study, we examine the potential gain of introducing the significance variable into the 
selection criteria for W — > ev analyses. The set of criteria employed is that of the recent mea- 
surement by the CMS collaboration of the W cross section 1251 , for which backgrounds were 
controlled using a stringent, 80% efficient, electron isolation criterion. Signal and background 
yields in that analysis were determined by a fit to the reconstructed Bj distribution, though 
because of the large backgrounds at small values of B t, the signal level is largely determined 
from the B t > 20 GeV region. 

One analysis option would be to relax the electron isolation from an 80% to a 95% efficient 



criterion and introduce B t significance to help control backgrounds. Figure 23 compares the 
efficiency for signal versus background in simulation for increasing minimum thresholds on: 
B t, with both the 80% and 95% electron isolation criteria applied; tSpF, again with both isolation 
criteria; or Bj/ \/YL Et, with the 95% isolation criteria. All efficiencies are measured relative to 
the signal or background yield obtained with the looser 95% electron isolation criterion ap- 
plied. (As a result, the tighter 80% criterion has an asymptotic value of approximately 84%.) 
Application of the tighter criterion changes the relative signal and background distributions 
for B t and B t significance compared to the looser criterion. When a minimum B t threshold 
is applied, the tighter isolation criterion provides a better signal to background ratio at low 
background levels than the looser criterion. Application of a minimum 5pf threshold with the 
looser criterion, however, outperforms all the other combinations for background rejection at a 
given signal efficiency. 

We note that in the calculation of the significance, the isolated signal electron candidate enters 
as an electron, and in particular with the resolution associated with an electron. This approach 
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was found to outperform the option where each event was treated as electron-free (as is the 
case for the dominant background). 



Figure 23 also shows that the 5pf distributions for W — > ev in data and simulation agree well. 
As expected, the backgrounds without genuine E j are compressed towards low values of 5pf 
while signal events having real tj extend to high values of 5pf- 
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Figure 23: (left) Efficiency curves for W — > ev signal versus backgrounds varying the minimum 
value of £x (solid lines), of Sp-p (dotted lines), and of Bt/ y/YL Er, (dot-dash line), with the 95% 
efficient (blue) or 80% efficient (red) electron isolation criterion applied, (right) Distributions 
for <Spf in candidate W — > ev events from data (points) and simulation (stacked histograms). 
The simulation components, from top to bottom, are signal (mustard) and backgrounds from 
jets (purple), 7+ jets (black), Z — > e + e~ (yellow), and W ± — > t ± v t (orange). The simulation is 
scaled by a fit to the data with floating normalizations for the signal and the total background. 



Figures 24 and 25 contrast the behaviour of signal and total background efficiencies for mini- 
mum £ x or 5pf thresholds for different numbers of interaction vertices (pile-up) in simulation. 
The jets and 7+ jets backgrounds, which have no genuine B7, dominate. The background con- 
tribution at higher Bj grows as pile-up increases, while the tSpp levels remain quite stable. As 
a result, a background subtraction based on extrapolation of E t will be sensitive to the mod- 
eling of pile-up, while one based on extrapolation of 5pf would not. As one can see from the 
signal versus background efficiency curves shown in Fig. 25 differentiation of signal from 



background degrades for both £x and 5pf as pile-up increases. Regardless of the amount of 
pile-up, however, <Spf always provides a superior signal to background ratio compared to tj. 
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Figure 24: Efficiency versus minimum threshold curves for W — > ev signal and for total back- 
ground for different numbers of interaction vertices with a minimum applied Bt threshold 
(left) and a minimum applied 5pf threshold (right). 
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Figure 25: Efficiency curves for W — > ev signal versus backgrounds varying the minimum 
value of B t (solid lines) and of 5pf (dotted lines) for events with two interaction vertices (blue) 
or at least five interaction vertices (red). 
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9 Conclusions 

In conclusion, we studied as measured by the CMS detector with pp collisions at a/s = 7 

— > 

TeV during the 2010 run of the LHC. We have determined the Bj scale, resolution, tails, and 

— * 

the degradation of the Bj performance due to pile-up. We find that we are able to remove most 
sources of anomalous energies that produce artificial, large Bj- The measured Bt scale agrees 
with the expectations of the detector simulation, but the resolution is degraded by 10% in data. 
CMS has three different algorithms for calculating £ T . Algorithms using tracker information 
have an improved resolution, and the use of a global particle-flow event reconstruction gives 
the best resolution. We find that pile-up interactions contribute to the degradation of the B t 
resolution but have little effect on the scale of a component of Bj- We also find that we can 
model its effects with a simple parameterization. 

One of the most important uses of B t is to distinguish between genuine B t produced by weakly 
interacting particles and artificial Bj from detector resolutions. An algorithm, called Bj signif- 
icance, for separating genuine B t from artificial B t, is shown to perform better than traditional 
variables such as Bj alone or Bj divided by the square root of the Ej. 

Because of the demonstrated good measurement of B t, the CMS detector is ready to be used 
for a variety of precision physics measurements, such as studies of the W boson, the top quark, 
and searches for new neutral, weakly interacting particles. 
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A Appendix: Optimization of £ T Corrections 

In this section, we describe the optimization of the parameters used in the type-I and type-II 
corrections to £t- 

For type-I £ j, the pj threshold used to select the jets that receive the correction was optimized 
by examining its effect on the £ t resolution and scale, and led to our choice of 10 GeV for PF £ j. 
We have optimized the pj threshold to obtain the best £ j scale and the best £ t resolution at the 
same time, under the constraint that very low-px jets should not be included in the calculation 
because their energy corrections have large uncertainties and can degrade the £ j performances. 
The use of low pj jets also makes the measurement more sensitive to the pile-up. 
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Figure 26: Response for PF £x versus 7 qj for various jet corrections thresholds for 7 candidate 
events with 1 primary vertex (left) and more than 1 vertex (right). 



Figure 26 shows the PF£j response versus qj for 7 candidate events for various jet thresh- 
olds for events with and without pile-up. As can be seen from these distributions, the type-I 
corrections substantially improve the £j scale, the dependence on the threshold is small for 
thresholds < 10 GeV, and is independent of the number of additional pile-up interactions. 



Figure 27 shows the ratio of the PF £ j resolution with type-I corrections to the PF £ j resolution 
without any type-I corrections for the parallel and perpendicular components of PF £ j versus 
qj for 7 candidate events for various jet thresholds and for events with and without pile-up. 
The optimal resolution for 1 vertex events is obtained for the chosen 10 GeV threshold, and for 
this threshold, the degradation with pile-up is not large. 

The same study has been performed in Z — > candidate events. The 10 GeV threshold is 

also the best threshold to combine the improvement of the £ j scale and £ j resolution in these 
events. 

Type-II corrections were not used for analyses done using the 2010 data. However, studies 
were done to optimize the parameters for future use. The type-II correction is obtained using a 
Z — > ee data sample. Z events constitute an ideal sample for determining the type-II correction, 
as Z's are generally produced with low qj and the recoil is often dominated by unclustered 
energy. We take the vector sum of momenta defined for all calorimeter towers or PF particles 
not corrected by the type-I correction as a single object denoted U. The measured value of U is 
defined as: 
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Figure 27: Ratio of PF B j resolution with type-I corrections to the the PF B j resolution without 
any type-I corrections for parallel (left) and perpendicular (right) components of PF B t versus 
7 cjj for various jet corrections thresholds for 7 candidate events with 1 primary vertex (top) 
and more than 1 vertex (bottom) after correction for scale, using scale corrections similar to 



those shown in Fig. 10 



^T,meas — B Tuncorr PT,jet,uncorr 2_j P T ,ele,meas/ (10) 

where B Tuncorr is the uncorrected Bt, PT,jets,uncorr is the momentum of uncorrected jets, and 
PT,eie,meas is the momentum of the measured electrons. The sum is over all jets with the corrected 
px > 20 GeV for Calo Bt and p T > 10 GeV for PF £ T . 

The type-II correction for the unclustered energy U meas is obtained by selecting events without 
any reconstructed jets using the correlation between U meas and the qj of the Z measured from 
the electrons. 

As the direction of the U meas m ay differ from the direction of qj due to noise, underlying 
event, etc., the parallel component of U meas projected on the direction of qj, Q meas n, is used 
for the derivation of the correction. The response of the unclustered energy is then defined 
as R(U r x,meas) = ^T,meas,|| /?T- The obtained correction factor, LT sca i e , is parametrized as 2.3 + 
1.9 exp(— 0.2 ■ lfr,meas) f° r CaloEj- For YYBj the obtained correction leads to a U sca \ e of 1.4. 
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These values obtained from data correspond well to the MC expectation. 

One of the ways to validate the type-II corrected £ j is to look at £ j decomposed into two com- 
ponents in dijet events. The decomposition is based on the dijet bisector axis, which divides 
the azimuthal opening angle between the leading two jets in two as illustrated in Fig. 28 Pos- 
itive J? tii points towards the smaller opening angle of the leading two jets, and positive £t± 
points towards the more central jet (smaller \t]\). The distribution of £t± is symmetric for all 




Figure 28: Kinematic definitions for dijet events. 

correction levels, while by definition there is a slight asymmetry in This is due to the fact 
that the bisector axis always points towards the opening angle of the leading two jets, while the 
£ t tends to point in the opposite direction. The type-I correction introduces a more significant 
asymmetry in the Bjn distribution, because it produces artificial tj in the direction opposite 
to the dijet opening angle. The type-II correction, however, calibrates the rest of the calorimeter 
energies, and makes the Etm distribution nearly symmetric again. This trend may be seen in 



Fig. 29 which shows the mean values of J?tii as a function of py VS = (pj' 1 + Pj t2 )/2 for uncor- 
rected^ type-I corrected, and type-II corrected Calo- and PF Bj in events containing at least two 
jets with p T > 40 GeV and \t]\ < 3. 
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Figure 29: Mean values of Bju as a function of average pj of the leading two jets (p T ) in dijet 
events compared with Monte Carlo simulation (dashed lines) for uncorrected, type-I corrected, 
and type-II corrected (left) Calo £ j and (right) PF £ j . The blue band on the distribution with 
type-I corrections represents the uncertainty due to uncertainties on the jet energy scale. The 
red band on the distribution with type-II corrections represents uncertainties due to uncertain- 
ties on the jet energy scale and statistical uncertainties due to the size of the Z — > e + e~ sample 
from which the correction for unclustered energies was derived. 
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